Anomalous Conductance Oscillations and Half-Metallicity in Atomic Ag-O Chains 
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Using spin density functional theory we study the electronic and magnetic properties of atomically thin, sus- 
pended chains containing silver and oxygen atoms in an alternating sequence. Chains longer than 4 atoms 
develop a half-metallic ground state implying fully spin polarized charge carriers. The conductances of the 
chains exhibit weak even-odd oscillations around an anomalously low value of 0.1 Go (Go = 2e^h) which co- 
incide with the averaged experimental conductance in the long chain limit. The unusual conductance properties 
are explained in terms of a resonating-chain model which takes the reflection probability and phase-shift of a 
single bulk-chain interface as the only input. The model also explains the conductance oscillations for other 
metallic chains. 
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PACS numbers: 73.40.Jn, 73.40.Gk, 73.63.Nm, 85.65.+h 

Atomically thin metallic chains can be formed by pulling 
a metal contact apart using e.g. a scanning tunneling mi- 
croscope or a mechanically controlled break junction 01 
Due to the low coordination number of the involved atoms, 
these ultimate one-dimensional quantum wires exhibit high 
mechanical stability and are chemically far more reactive than 
their bulk counterparts |3, 4, 5]. This, combined with their 
ability to sustain extremely large current densities, makes 
atomic chains interesting from both a fundamental and a tech- 
nological point of view. 

At sufficiently low temperatures the conduction electrons 
in a metal can travel hundreds of nanometers without losing 
their quantum mechanical phase 111]. Under such conditions, 
the wave nature of electrons becomes important and manifests 
itself directly in macroscopic quantities such as the electrical 
conductivity. 

As first predicted by density functional theory (DFT) calcu- 
lations, the conductances of suspended atomic chains of some 
elements like C, Na, and Au oscillate with a period of two 
as the number of atoms, N, in the chain is varied |0, 
This behavior can be understood from a local charge neu- 
trality condition which implies that the Fermi level must be 
aligned with the center of a resonance for odd N and lie 
mid between two resonances for even N, respectively. De- 
viations from these even-odd oscillations naturally occur for 
metals with more complex valence configurations such as Al 
and Pt |8, 9]. Common to the oscillations observed for all the 
homogeneous metal chains is that the conductance maximum 
(per contributing channel) is always close to IGq while the os- 
cillation amplitude varies significantly from ~ O.OSGq in the 
case of Au |9, 10[] to ^ O.SGq in the case of Al and C chains 
. We show here, that this behavior can be understood in 
terms of the reflection probability and phase shift at a single 
bulk-chain interface. 

In a recent experiment Thijssen et al. found that the pres- 
ence of oxygen greatly enhances chain formation when a sil- 
ver contact is broken at low temperature [12]. It was pro- 
posed that the resulting chains, which can exceed 2 nm in 
length, consist of Ag and O atoms repeated in an alternating 



fashion [131], consistent with previous DFT calculations for 
Au-O chains ^ [sf] • Interestingly, the conductance averaged 
over many chains and plotted as a function of chain length 
showed an exponential-like decrease from IGq to O.IGq af- 
ter which the conductance stayed constant over almost 1 nm. 
This is indeed a surprising result: the initial decrease in con- 
ductance suggests a chain with a band gap, however, this is not 
in line with the saturation at a finite conductance for longer 
chains. On the other hand, a length-independent conductance 
of 0.1 Go does not comply with the results for metallic chains 
which exhibit conductance oscillations with a maximum on 
the order of IGq. We mention that silver, gold, and nickel 
contacts with one or two oxygen atoms adsorbed have previ- 
ously been studied. ]j5i [13. US 13] In the following we 
show that the conductance behavior follows from the simple 
phase-shift model. 
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FIG. 1: (Color online.) Calculated conductance as a function of chain 
length (blue squares). The length of the chain is given in units of 
3.8A corresponding to the length of a Ag-0 unit. Also shown is 
the result of the resonating-chain model (red triangles) and a fit to 
the experimental data of Ref. 12. The supercell used to model the 
suspended chains is also shown. 

In this Letter we present spin DFT calculations showing 



2 



that alternating silver-oxygen chains (ASOCs) exhibit small- 
amplitude (< O.OTGq) conductance oscillations around a 
value of O.IGq for long chains. The origin of these anoma- 
lous conductance oscillations is traced to an additional phase 
picked up by an electron being reflected at the end of the 
chain. Furthermore, we find the ASOC to be half-metallic 
meaning that only electrons with a given spin direction can 
pass through the chain making it a perfect spin-valve. 

To describe the chains we use the supercell shown in the 
inset of Fig. [T] It contains the oxygen terminated ASOC 
suspended between four-atom silver pyramids, which are at- 
tached to (1 1 1) silver surfaces. Periodic boundary conditions 
are imposed in all directions. Both the silver pyramid and the 
silver-oxygen chain have been relaxed using the Dacapo plane 
wave DFT code ifisl] to obtain the most stable geometry ifioll . 

Assuming phase-coherent transport, the conductance of the 
chains in the limit of low-temperature and -bias, is given by 
the Landauer-Biittiker formula, G — GoT{eF), where T(e) is 
the energy-dependent transmission function. The latter is cal- 
culated using the Non-equihbrium Green's function (NEGF) 



method described in Refs. 11201 12111 where the Kohn-Sham 
Hamiltonian of the relaxed chains is obtained from the Siesta 
code 122(1 . 

In Fig.[T]we show the calculated conductance (blue squares) 
of the relaxed oxygen-terminated ASOCs as a function of 
the number of oxygen atoms, N, in the chain. We note in 
passing that the calculated conductances of silver-terminated 
chains are significantly lower than the experimental value of 

0. 1Go [13], but these chains turn out to be less stable than 
oxygen-terminated chains and will probably break in the Ag- 
Ag bond when elongated. As indicated on Fig. [T| we have 
considered relaxed chains with a length up to 23.0 A, with the 
Ag-O bond length lying in the range 1.95 — 2.10 A. For com- 
parison the bond length of the infinite linear alternating Ag-0 
wire is found to be 1.95 A. 

Fig. [T also contains a fit to the experimental data presented 
in Ref. | l2l. The experimental data is obtained by averag- 
ing the conductance trace of thousands of chains, and fea- 
tures such as small-amplitude conductance oscillations are not 
visible. For long chains, however, there is excellent agree- 
ment between the measured and calculated (average) conduc- 
tance. The significantly larger conductances found in the ex- 
periments for shorter chains are presumably due to contribu- 
tions from pure Ag chains which have a conductance around 
IGo. 

The calculated conductance of a single oxygen atom, N — 

1, is 0.3Go which is in good agreement with the calculated 
results in Ref. [14]. For > 1 the conductance starts to 
oscillate with a period of two Ag-O units around an average 
value of 0.1 Go. For = 1 we find a non-magnetic ground 
state, while for > 1 the ferromagnetic state is energeti- 
cally favored. For large N the electronic structure of the chain 
converges toward that of the infinite alternating Ag-O chain 
which has a ferromagnetic ground state with a magnetic mo- 
ment of l/is per Ag-O unit mainly localized at oxygen. We 
find the energy gain with respect to the non-magnetic state to 



be 0.12 eV per Ag-O unit. This breaking of the spin symmetry 
is the reason for the lower conductance found for the chains 
with N > I. 



(a) 



2 
1 
2 

c 1 

o 

05 2 
05 

E 1 

05 

§ 2 





/ \ -- Majority 
/ V — Minority 


'1 

'» M r 




A 




A, A 




Vi 


A A-- 







3 m 



o 

CO 

03 
X 

o 



05 



-0.8 -0.4 0.4 

e-£p(eV) 



0.8 



(b) 



> 

CD 



OJ 

I 



-2 



1 1 1 
: -- Majority 


1 , 1 , - 


- — Minority 


in=±1- 




^^^^^ " " 




y"^ . ^"^ " " 







0.1 



0.2 0.3 

k {2n/a) 



0.4 



0.5 



FIG. 2: (Color online) (a) Transmission functions for the majority 
spin (dashed line) and minority spin (full line) for Ag-O chains con- 
taining 2-6 oxygen atoms, (b) Band structure of the infinite Ag-O 
chain. The band crossing the Fermi level is doubly degenerate and 
the Bloch states have m — ±1 symmetry with respect to rotations 
around the chain axis. 



To gain more insight into the nature of the conductance 
oscillations, we show in Fig. |2^ how the energy dependent 
transmission functions change with N. The number of oxy- 
gen atoms in the ASOC varies from iV = 2 in the top panel 
to iV = 6 in the bottom panel. We first notice that only mi- 
nority spin states are present at the Fermi level which means 
that the current will be fully spin polarized. The peaks in the 
transmission function coincide with resonances of the chain 
which are broadened by coupling to the contacts. The num- 
ber of resonances increase linearly with N and their width 
decrease as 1/iV due to a reduction in the overlap between 
the levels on the chain and the states in the contact. The res- 
onances have the same character as the Bloch states forming 
the valence band of the infinite Ag-O chain whose band struc- 
ture is shown in Fig. The spin polarized valence band is 
two times degenerate with angular momentum m — ±1 and 
Ag(4d)-0(2p) character. 



Clearly, the conductance oscillations arise because the 
Fermi level intersects the nearest resonance closer to its cen- 
ter for even N than odd N. In contrast to the situation for 
homogeneous chains, however, the transport is never "on res- 
onance" but always takes place via a resonance tail. In par- 
ticular this cannot be explained by charge neutrality. Indeed, 
from Fig. it can be seen that the resonances are always 
almost completely empty or filled, and since each resonance 
can accommodate two electrons (due to orbital degeneracy), 
it seems that all chains accommodate an even number of elec- 
trons. But each Ag-O unit contains 5 valence electrons, and 
local charge neutrality would therefore imply a half-filled res- 
onance for every second Ag-O unit added to the chain. On the 
other hand, the small-amplitude oscillations around the exper- 
imental saturation value of O.IGq arise exactly because the 
Fermi level always intersect the tail of a resonance. 

In order to rationalize the conductance behavior of the 
ASOCs, as well as that of other chain systems, we propose the 
following simple resonating-chain model, which is exact for 
coherent transport above a certain chain length. An electron 
approaching the chain from the left lead will have a certain 
transmission probability Ti for entering into the chain. Inside 
the chain the electron can propagate back and forth between 
the two (identical) contacts, so that every time the electron 
impinges on one of the contacts it is reflected with a proba- 
bility i?i = 1 — Ti and it furthermore picks up a phase shift 
which we denote by 01. We furthermore assume that the elec- 
tron in the infinite chain is characterized by a single Fermi 
Block wave vector k. The interference between all the re- 
flected waves then leads to a total transmission probability T 
through the chain system of 



l + 4(ii§7)7sin2(fcL + </>!)' 

where L denotes the distance between the two contacts. This 
expression is exact for coherent transmission in the limit 
where the reflections at the two ends can be considered in- 
dependent, i.e. it is required that the potential in the central 
part of the chain is unaffected by the two contacts. It should 
be noted that there is a certain arbitrariness in the definition 
of the distance L between the two contacts, however, different 
choices for L also gives different results for the phase shift (pi 
so that Eq.[T]still holds. We take the points of reflections to be 
at the tip atoms of the contacts. 

For a given reflection probability Ri the total transmission 
is seen to vary with the chain length between a maximum of 1 
and a minimum of((l — i?i)/(l + i?i))^. However, due to the 
phase shift and the discrete nature of the length of the chain 
the transmission usually oscillates over a more narrow region. 

We have determined the parameters -Ri and cf>i for the 
ASOCs as well as for chains of Al and Au. The param- 
eters are obtained by considering the reflection of an elec- 
tron in a semi-infinite chain impinging on a contact (see 
Fig.O. We determine the scattering state as described 
in Ref. ll24ll . and the phase shift is then determined by pro- 
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FIG. 3: (Color online) Isosurface of the absolute value of the eigen- 
channel state vfs for a Ag-O chain. The state is a superposition of 
incoming and reflected chain Bloch states. 



jection onto the incoming (k) and outgoing (— fc) states as 

(f>i =arg((?A_fc|«',)/(V^fc|'J',)). 



System 


Ri 




kL 


Tfnax 


^Tosc 


Ag-O/Ag 


0.64 


-0.2l7r 


Ntv/2 


0.12 


0.05 


A1/A1(111) 


0.32 


0.927r 


(iV + l)7r/4 


0.85 


0.57 


Au/Au(100) 


0.004 


-0.367r 


{N + 1)tt/2 


0.997 


0.010 



TABLE I: Calculated chain-contact reflection parameters for a num- 
ber of chain systems. 

In Table U the calculated reflection parameters for the 
ASOCs and for mono-atomic Al and Au chains are shown to- 
gether with the obtained maximal transmissions and the mag- 
nitudes of the oscillations. In all three cases the model repro- 
duces the results of full DFT calculations for varying chain 
lengths. 

The model results for the ASOCs are shown in Fig.[T]to co- 
incide with the full calculations for chains containing two or 
more oxygen atoms. It should be noted that the low conduc- 
tance of about O.lGo comes about not so much because of 
the reduced reflection probability (_Ri = 0.64) which would 
still allow for a total transmission of 1 if a resonance condi- 
tion could be met in the chain. Rather it is the phase shift of 
01 = — 0.2l7r K, — 7r/4 which leads to destructive interfer- 
ence within the chain. The actual size of this phase shift de- 
pends on both the electronic structure of the chain and of the 
contacts and must be obtained from a full calculation. How- 
ever, we note that a phase shift of 7r/4 corresponds to a length 
of half a unit cell in the chain or equivalently an oxygen-silver 
distance, and this is in good accordance with Fig.[3]where the 
scattering state is seen to carry weight on the Ag tip atom and 
therefore is rather reflected off the three-atom Ag layer below. 
We also mention that Eq.[T] reproduces the energy-dependent 
transmission (Fig.|2^) quite well including the resonant struc- 
ture, thus implicitly resolving the charge neutrality paradox - 
odd numbers of chain electrons still doubly-occupy an integer 
number of resonances when the effects of the phase shift at 
the boundaries are included. 

The results of the resonating-chain model for the conduc- 
tances of Al and Au chains are shown in Fig. |4] together with 
the results of earlier full DFT-transport calculations. The 
agreement is very good even for rather short chains. In the 
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FIG. 4: (Color online) (a) Conductance for Al as function of the 
number of chain atoms, A'^. The first-principles result (squares) is 
taken from Ref. 8 while model result (triangles) is obtained from 
Eq. [T] with parameters derived using the same procedure as for the 
Ag-O system, (b) Conductance length dependence for Au with the 
first-principles result taken from Ref.fiol 



case of Al the phase shift is close to tt in excellent agree- 
ment with the resonant-level model proposed in Ref. |@]. The 
resonant-level model takes as a starting point the isolated 
chain (of length N) which is then coupled weakly to the leads, 
and this boundary condition corresponds exactly to a phase 
shift of TT. The small deviation of the phase shift from tt ac- 
counts nicely for the fact that the conductance does not peak 
at 26*0 (the factor of two coming from two available channels) 
but only at I.TGq. 

For Au we find the phase and amplitude of the even-odd 
conductance oscillations in agreement with the first-principles 
calculations in Ref. 10. We note, that the phase of the conduc- 
tance oscillation is not explained by either a charge neutral- 
ity argument or a resonant model [8]. However, the phase of 
the oscillation can be related to our calculated phase shift of 
— O.SGtt, while the small oscillation amplitude can be traced 
to an almost perfectly transparent chain-bulk interface with a 
reflection coefficient of 0.004. 

In summary, we have presented first principles conductance 
calculations showing that alternating Ag-O chains suspended 
between silver bulk contacts are half-metallic and can have 
an average conductance of O.IGq as found in a recent experi- 
ment llZ]. In fact, the conductance oscillates with a small am- 
plitude and a period of two Ag-O units as the chain length is 
varied. The oscillations can be understood from a resonating- 
chain model, and are fully characterized by only two parame- 
ters: The reflection probability, and reflection phase-shift, 
01, of a single bulk-chain interface. By extracting these two 
parameters from the DFT calculation, quantitative agreement 
between the full calculations and the model is obtained. 
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